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1. Introduction 


1.1 Background Information 


In computational complexity theory, a problem is assigned to the NP (non-deterministic 
polynomial) class if it can be verified in polynomial time. The Travelling Salesman 
Problem (TSP) is potentially the most famous optimisation problem and it falls under the 
NP-hard category since the existence of a polynomial-time solution for it implies the 
existence of a polynomial-time solution for every problem in NP.! The TSP consists of 
determining the shortest tour to complete a Hamiltonian Cycle — a path through a graph 
that starts and ends at the same vertex, including every other vertex exactly once; an 


example is shown in Figure 1.2 


Figure 1: A Hamiltonian cycle with 8 vertices? 


'NP-hard problems and approximation algorithms - University of Texas at... Available at: 
https://personal.utdallas.edu/~dxd056000/cs6363/unit5.pdf (Accessed: 17 May 2023). 

2 Black, P.E. (2020) Hamiltonian cycle, Dictionary of Algorithms and Data Structures Available at: 
https://xlinux.nist.gov/dads/HTML/hamiltonianCycle.html (Accessed: 17 May 2023). 

3 Chatting, M. (2018) 'A Comparison of Exact and Heuristic Algorithms to Solve the Travelling Salesman Problem', The 

Plymouth Student Scientist, 11(2), p. 53-91. 


1.2 Contextual Significance 


The TSP has stimulated the development of various problem-solving techniques, 
algorithms and innovative mathematical models that can be applied beyond its immediate 
problem space. For instance, this can be applied to the sphere of network and hardware 
optimisation where the most efficient route for data transmission ought to be found. The 
formulation as a TSP essentially provides the simplest way to solve problems arising in 
many different contexts, including computer wiring, vehicle routing, clustering, and job- 


shop scheduling. 
1.3 Scope of Research 


There are two variations of the TSP: asymmetric Travelling Salesman Problem (ATSP), 
where the distance from node A to B differs to that from B to A, and symmetric Travelling 
Salesman Problem (STSP), where the graph is undirected.^ Hence, regarding the STSP, 
it can be said that cj; = cj; this simply states that regardless of the direction of travel, the 
cost (or distance) between city i and city j is constant.” Numerous algorithms have been 
generated to approximate a solution to the TSP in a feasible time span since finding the 


optimal solution for large problem instances is computationally challenging. 


The aim of the paper is to determine the most efficient algorithm for solving the STSP by 
analysing both the time complexity and accuracy of the algorithms. In this case, the most 
efficient algorithm can be quantified as the one that has the greatest accuracy and 


shortest execution time. 


^ Deep, Kusum., & Mebrahtu, Hadush. (2012), "Variant of partially mapped crossover for the Travelling Salesman 
problems." International Journal of Combinatorial Optimization Problems and Informatics, Vol.3, num.1, pp.47-69. 
ISSN: 2007-1558 

5 Ibid. 


1.4 Experimental Overview 


This paper specifically focuses on comparing the Elitist Genetic Algorithm to Ant Colony 
System, evaluating which of the two is most efficient at solving the STSP. Random data 
sets of increasing size will be used to collect a set of execution time periods for each of 
the algorithms. Furthermore, the accuracy of the algorithms will be obtained by comparing 
the shortest distance calculated to that determined by a control algorithm; for this 
experiment, the Brute Force Algorithm will be used - an exact algorithm, so the optimal 
solution is guaranteed to be found, with time complexity of O(n!).* These two factors were 
then analysed to answer the question: "How does the Elitist Genetic Algorithm 
compare to Ant Colony System in terms of time complexity and accuracy when 


attempting to solve the Travelling Salesman Problem?" 


6 Chase, C.et al.(no date)An Evaluation of the Traveling Salesman Problem. Available at: 
https://scholarworks.calstate.edu/downloads/xg94hr81g#:-:text=. (Accessed: 14 July 2023). 


2. Research 


2.1 Genetic Algorithms 


Genetic Algorithms (GA) are adaptive, stochastic, metaheuristic search algorithms that 
are a subclass of evolutionary computing used to solve combinatorial optimisation 
problems.” Combinatorial optimisation is the process of solving for the optimal solution of 
a finite data set by using combinatorial techniques.? Metaheuristics are algorithmic 
concepts that define heuristic methods applicable to a number of different problems.? 
Darwin's theory of evolution (survival of the fittest and natural selection) forms the 


backbone of GAs; Figure 2 describes the steps of this process. 


Fitness calculation 


Output result 


Crossover 
Repair and mutation 
Preserve elite 


Figure 2: Flowchart for the Elitist Genetic Algorithm” 


7 Genetic algorithms (2017) Scribd. Available at: https://www.scribd.com/document/351623322/Genetic-Algorithms# 
(Accessed: 15 June 2023). 

8 Ibid. 

9 Dorigo, M. and Stützle, T. (2004) ‘The Ant Colony Optimization Metaheuristic’, in Ant colony optimization. 
Cambridge, MA: MIT Press, pp. 25-26. 

'0 Singh, V.K. and Sharma, V. (2014) ‘Elitist genetic algorithm based energy balanced routing strategy to prolong 
lifetime of wireless sensor networks', Chinese Journal of Engineering, 2014, pp. 1—6. doi:10.1155/2014/437625. 


6 


These algorithms are more robust and can navigate through larger data sets whilst finding 
an optimal solution within a sensible timespan. Chromosomes are used to represent n 
number of genes (cities) in the order in which they are visited, typically as a binary string. 
Figure 3 shows a chromosome of a 5-city tour, where the salesman begins at city 3, then 


travels to city 4 and so on. 


3.2|4. 1,5 
Figure 3: Chromosome of a 5-city tour 


2.1.1 Population Initialisation 


In terms of the TSP, the city tour is referred to as the population. The first process of a 
GA is to initialise the population — the cities can be randomly generated or set before. If 
the initial population size is too small, diversity is prohibited which forces some 
optimisation routines to converge too quickly, causing the population to become 
homogenous.'' However, if the population size is too big, it would take a long time for the 


optimisation routine to converge. '2 


2.1.2 Fitness Evaluation 


As each city is situated in a 2D plane, their position can be given by (xj, yj). Using 


Equation (1) and the position coordinates of two cities, the cost c;;, in other words the 


j: 


distance between the ith and jth city, can be found: 


ej = Ju — x)? + (y, — yy? 8 (1) 


11 Morris, A.T. (1998) Optimization of the Traveling Salesman Problem and Multivariate Real-Valued Functions using 
a Genetic Algorithm. dissertation. 

12 Ibid. 

13 Deep, Kusum., & Mebrahtu, Hadush. (2012), op. cit. 


The distances between the cities can be represented in an n x n cost matrix. Below is an 
example of a symmetric 5-city tour. 
City 1 City 2 City 3 City 4 City 5 

Cityl r 0 8 14 23 2 

Cty2| 8 0 26 2 19 

Cty3 114 26 0 2 23 

Ciy4 |23 2 2 0 11 

Cty5 L2 19 23 11 0 
Each chromosome is then evaluated using a fitness function and is assigned a fitness 
value. The fitness value is determined based on the distance between the two cities using 
Equation (2): 


Fitness = Yj-., distance(x; to xj44) !^, (2) 


where n is the number of cities and i is the city index value. 


The shortest distance is given the highest fitness value. 


2.1.3 Selection and Elitism 


A random sample of chromosomes is placed in a mating pool via the reproductive 
process, selection, by which the fitter chromosomes, those with shorter city tours, are 
more likely to be reproduced to the next generation."® Elitism is also considered part of 
the selection process for this paper, meaning it is guaranteed that the best chromosome, 


the one with the shortest distance, is copied to the next generation.'® 


14 Deep, Kusum., & Mebrahtu, Hadush. (2012), op. cit. 

15 Morris, A.T. (1998) Optimization of the Traveling Salesman Problem and Multivariate Real-Valued Functions using 
a Genetic Algorithm. dissertation. 

16 Ibid. 


2.1.4 Evolution through PMX Crossover 


Crossover is a genetic operator by which design characteristics between two parent 
chromosomes, chosen randomly from the mating pool, are exchanged to form two new 
superior offspring.*” For this experiment, the Partially Mapped Crossover (PMX) operator 
will be used as one of the most effective and popular. It is a modification of the basic two- 
point crossover, however, uses an additional mapping relationship to avoid duplicate 
values in the offspring that often lead to infeasible results.'® PMX falls into the category 
of Inventing Specialised Operators - meaning only valid chromosomes are generated 


(cities are not missing or repeated).'? 


Below is an example PMX crossover, where P, and P, are random, parent chromosomes 


of an 8-city tour. Offspring O, and O, are formed which each represent a new city tour. 


P, = (41257368) 
P, = (15836247) 
A substring is selected using two random crossover points (marked with "|"): 
P, 2(412|573]68), 
P, 2(158|362]|47). 
A Two-Point Crossover is performed: 
0, = (xxx|362|xx), 


0, = (xxx|573|xx). 


17 Hasançebi, O. and Erbatur, F. (2000) ‘Evaluation of crossover techniques in genetic algorithm based optimum 
structural design’, Computers &amp; Structures, 78(1—3), pp. 435—448. doi:10.1016/s0045-7949(00)00089-4. 

18 Deep, Kusum., & Mebrahtu, Hadush. (2012), op. cit. 

19 Üçoluk, G. (2002) ‘Genetic algorithm solution of the TSP avoiding special crossover and mutation’, Intelligent 
Automation &amp; Soft Computing, 8(3), pp. 265-272. doi:10.1080/10798587.2000.10642829. 
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The mapping systems are determined: 
53,706,392. 
Bits that are not conflicting are filled: 
0, = (41x|362]|x8), 
0, 2(1x8/|573]|4x). 
Using the mapping relationships, the offspring can be fully filled: 
0, = (415|362|78), 


0, 2(128|573|46). 


2.1.5 Evolution through Mutation 


Mutation is a unary genetic operator which produces spontaneous, random changes on 
one parent chromosome.”° Only one type of mutation, the swap mutation, will be used to 
ensure reliability in the results as this acts as a control variable. Two genes (cities) are 


selected at random, and their positions are swapped. 


An example swap mutation is shown, where P, is the parent chromosome of an 8-city 
tour. 
P, = (41257368) 
Two cities are chosen: 
P, = (41257368). 
Their values are interchanged: 


0, = (43257168). 


20 GeeksforGeeks. (2018). Mutation Algorithms for String Manipulation (GA). [online] Available at: 
https://www.geeksforgeeks.org/mutation-algorithms-for-string-manipulation-ga/. 
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2.2 Ant Colony Optimisation 


Ant Colony Optimisation (ACO) is another stochastic, population-based, metaheuristic 
search algorithm that simulates the foraging behaviour of ants to solve combinatorial 
optimisation problems.? The concept of using swarm intelligence, the collective 
behaviour of decentralised, self-organised natural or artificial systems, was first 


introduced by Gerado Beni and Jing Wang in 1989.7? 


2.2.1 Real Ant Behaviour 


Ants use stigmergy meaning they indirectly communicate with each other by altering their 
surrounding environment. They have collective intelligence as they lay a pheromone trail 
while searching for food to communicate with each other to find the shortest path; other 
ants can sense this chemical, influencing their choice of path. Pheromone is a particularly 
volatile substance that starts to evaporate after the ant marches over the path. Hence, for 
shorter paths, the pheromone density remains high as pheromone accumulation is 


faster.?? Figure 4 shows how the distribution of pheromones is dependent on path length. 


A 


Nest Food 


[| [| 
Figure 4: Pheromone trails of path A and B^ 


21 Ahmed, Z.E.et al. (2020) ‘Energy optimization in low-power wide area networks by using heuristic 
techniques', LPWAN Technologies for loT and M2M Applications, pp. 199—223. doi:10.1016/b978-0-12-818880- 
4.0001 1-9. 

22 Ibid. 

23 Ranjith, K.A. (2010) Ant Colony Optimization. rep., pp. 1—16. 

24 Nguyen, K.-H. and Ock, C.-Y. (2011) ‘Word sense disambiguation as a traveling salesman problem’, Artificial 
Intelligence Review, 40(4), pp. 405-427. doi:10.1007/s10462-01 1-9288-9. 
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Initially, there is an equal probability p that an ant will travel via path A or B when in search 
of food. As path A is shorter than path B, in a specific time period t, path A will be travelled 
more times; therefore, path A will have a higher pheromone density. As t increases, more 
ants will follow path A whilst the pheromone trails in path B will all evaporate - eventually 


all ants follow path A.?® 


2.2.3 Ant Colony System 


There are many variations of ACO algorithms including Rank-Based Ant System 
(ASrank), Max-Min Ant System (MMAS) and Ant Colony System (ACS). In this paper, the 
ACS algorithm will be investigated - a set of cooperating agents, ants, indirectly 
communicate with each other through the deposited pheromones on the edges of the 
TSP graph whilst finding the optimal solution.? All ants perform the local pheromone 
update after every step rather than after a completed tour meaning the next edge is 
chosen purely based on the updated pheromone value. The process stops when the best 
solution is found or there are no more pheromone updates.” During ACS, an ant 
completes a tour around the map n number of times. After every iteration, the ant's global 


memory is reset; meaning, they have no knowledge of the journey they took prior. 


25 Ranjith, K.A. (2010), Op. cit, 

26 Dorigo, M. and Gambardella, L.M. (1997) “Ant Colony System: A cooperative learning approach to the traveling 
salesman problem’, IEEE Transactions on Evolutionary Computation, 1(1), pp. 53-66. 
doi:10.1109/4235.585892. 

27 Mulani, M. and Desai, V.L. (2018) 'Design and Implementation Issues in Ant Colony Optimization', in /nternational 
Journal of Applied Engineering Research. 16th edn. Research India Publications, pp. 12877-12882. 
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Figure 5 below describes the steps of an ACS where the evaluation stop condition is when 


n number of iterations have taken place. 


Initial Population 


| Randomly Ant | 
Located 


EN 


Each Ant 4 
< Randomly Chose a >< 
path 


A Y 
[vee Pheremone | | search process 
| tria J continue 


| Ph = 
BE 
Y 


Evaluation 
Stop y 
Condition 


bo 
| Results | 


Figure 5: Flowchart for ACS?? 


2.2.3 ACO Algorithms for the TSP 


When solving the TSP, this algorithm assumes that ants are always able to determine the 
shortest path to the food sources by detecting the pheromones laid by other ants. In other 
words, cities further away are less visible meaning there is a lower probability of being 
chosen. The greater the intensity of the pheromone trail, the greater the probability that 


the ant will choose that edge. 


28 M. Almufti, S., Boya Marqas, R. and Ashqi Saeed, V. (2019) ‘Taxonomy of bio-inspired optimization 
algorithms’, Journal of Advanced Computer Science &amp; Technology, 82), p. 23. 
doi:10.14419/jacst.v8i2.29402. 
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Pheromone trails describe the desirability of an ant visiting node j after node i. As evident 
in Equation (3), the heuristic desirability, 7;;, of an ant going from node i to node j is 


inversely proportional to the distance, d;;, between the nodes: 


j> 


nij = 1/d;;.? (3) 


At each node, the ant plans based on its local memory — this stores information about the 
adjacent nodes. Nodes that have not been visited by the nth ant are defined as allowed, . 


The probability of the nth ant choosing one node is given by Equation (4): 


[ij (t)]* [nig l? : 
= eae JE allowed 
pij (t) = Eneallowedpltij(t)l1 [nj] J allowed, (4) 
0 else, 


where 7;; is the evaporation rate, determining the amount of pheromone present between 
node i and j, a = 0 is a parameter in control of the influence of 7;;, nij is the desirability 
of the transition from i to j and f > 1 is a parameter in control of the influence of r;;.?' 
These parameters must be selected properly otherwise the convergence speed may be 
too high causing the algorithm to fall rapidly into local optima. Likewise, if the parameters 


are set so that the convergence speed is slow, time complexity increases. 


30 Dorigo, M. and Stützle, T. (2004) 'Ant Colony Optimization Algorithms for the Traveling Salesman Problem', in Ant 
colony optimization. Cambridge, MA: MIT Press, pp. 67-68. 

31 Danu, M.S. (2013) Ant colony optimization algorithms, Scribd. Available at: 
https:/Awww.scribd.com/document/136679005/Ant-colony-optimization-algorithms# (Accessed: 27 June 2023). 
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2.3 Parameters for Analysis 


2.3.1 Time Complexity 


The first parameter analysed for each algorithm is the time complexity. For a given 
problem, the time complexity is defined as the maximum time the algorithm requires to 
find a solution for each possible input size, n.?? This is alternatively referred to as the 
worst-case time complexity. Big-O notation is typically used to describe time complexity 


for the function f (n) as it gives asymptotic upper bounds for the worst-case scenario: 


fm) =0(9m)) forn>® and f(n),g(n) ER 


where g(n) represents the big-O notation. Equation 4 is only valid if there exist 


constants c and ny such that 


FOIS clg(| for all n> ng. 


This effectively means that the big-O notation, denoted by g (n), is always greater than or 


equal to the number of steps. Determining the time complexities of different algorithms 


enables their efficiency to be compared and analysed. 


*? Big O notation - mit - massachusetts institute of technology (no date) Big O Notation. Available at: 
https://web.mit.edu/16.070/www/lecture/big o.pdf (Accessed: 26 June 2023). 

33 Ibid. 

34 Ibid. 
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For example, as seen in Figure 6, an algorithm with a time complexity of O(logn) is 


significantly more efficient than one which has a time complexity of O(n!). 


a 
< 
2 
x 
S 
à 
o 


Elements 
Figure 6: A graph to show different time complexities? 


The execution time of NP-hard problems increases exponentially with the data size; 
consequently, so does the time complexity. For instance, the TSP has a time complexity 
of 0 (n!). Although finding the solution for low values of n is manageable, this is not the 
case for most scenarios. Hence heuristic approaches are used which have a better time 


complexity even though the solution may not be optimal. 


2.3.2 Accuracy 


Accuracy is the second parameter investigated in this paper as it is used to validate and 
assess the performance of a specific algorithm. It is defined as an evaluation metric that 


measures the closeness of the obtained values to the accepted or correct value. For this 


35 Prado, K.S. do (2020) Understanding time complexity with python examples, Medium. Available at: 
https://towardsdatascience.com/understanding-time-complexity-with-python-examples-2bda6e8158a7 
(Accessed: 26 June 2023). 
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experiment, the accuracy states the error between the optimal path generated by an 
algorithm and the actual, shortest path available. As mentioned previously, exact 
algorithms, including the Brute Force Algorithm, always find the shortest route for the TSP 
despite being extremely inefficient. Hence, the percentage accuracy of the Elitist GA and 


ACS will be found using Equation (5). 


shortest path generated by control brute force algorithm 


percentage accuracy = x100 (5) 


shortest path generated by the elitist GA or ACS 


2.4 Hypothesis 


Although the two algorithms fall under the same combinatorial optimisation category (as 
they are both population-based, metaheuristic and stochastic), this paper hypothesises 
that the elitist GA will have a lower time complexity compared to ACS. This is because 
the additional parameter, elitism, ensures that the best chromosomes are always 
transferred to the next generation, so the optimal solution should be obtained faster. 
However, ACS should be more accurate since there are several parameters that can be 
adjusted to acquire the desired convergence speed — this includes the desirability and the 
constants a and $. Therefore, the hypothesis for this investigation is: When using the 
elitist GA and ACS to find the shortest tour for the TSP, the Elitist GA should have 
a lower time complexity whilst the solution generated by ACS should have a higher 


level of accuracy. 
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3. Experimentation 


3.1 Methodology 


3.1.1 Variables 


The two independent variables for this investigation are the algorithms used, Elitist GA or 
ACS, and the data set size. As evident in Appendix 4, the number of cities to visit 
increases by 4 each time up until 20 as this provides a variety of results that can later be 


interpreted. 


As stated in the research question both the accuracy and efficiency of each algorithm will 
be measured. Hence, the two dependent variables are the shortest distance calculated 
and the time taken, in seconds. This enables the accuracy of each algorithm to be 


calculated using Equation (5) as well as the time complexities to be compared. 


One key control variable is the type of control algorithm used since it will produce the 
shortest distance of each city tour which will then be compared to those calculated by 
ACS and GA . This way the accuracy can be calculated. Hence, it was decided that the 


Brute Force algorithm will act as the control and the Python code to run this will be 


obtained from https://github.com/Joseph bakulikira/Traveling-Salesman-Algorithm. 
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Table 1 discusses all the control variables for this experiment. 


Controlled Variable 


Why is it controlled? 


How is it controlled? 


Computer and 
Operating System 


The hardware can impact the 
speed of each algorithm, 
depending on factors such as 
processor speed and memory. 
Therefore, these elements 
must be kept constant to 
ensure repeatability. 


The same hardware was 

used throughout the 

experiment: 

- Computer model: 21.5- 
inch 2019 iMac 

- Processor: 3GHz 6-Core 
Intel Core i5 

- Memory: 16GB 2667 
MHz DDR4 

- macOS: Ventura 13.4.1 


All applications will be closed 
when running code to keep 
the amount of RAM available 
constant. 


Integrated 
Development 
Environment (IDE) 


The specifications and differing 
features of an IDE can have an 
effect on the results obtained. 


The same IDE was used: 
- Type: Visual Studio Code 
- Version: 1.80.1 


Data points for each 
data set 


If different points are used, the 
shortest distance will be 
different so the accuracy of 
each algorithm cannot be 
calculated. 


The points used were 
chosen in the preliminary 
experiment and used 
throughout the experiment 
(see Appendix 4) 


Parameters for each 
algorithm 


The specific parameters of 
each algorithm have effects on 
factors like termination criteria, 
convergence speed and the 
exploration of space. 


The mutation rate for the 
Elitist GA was 0.0196. 

The parameters for ACS 
were fixed: 

- ait 

n8 


Initial population size 


The population size influences 
the convergence speed and so 
must be controlled. 


The population size for both 
algorithms was set at 150. 


Python code for ACS, 
GA, and Brute Force 


The same code must be used 
for each algorithm to ensure 
reproducibility and consistency 
in the results. 


All algorithms were taken 
from: 
https://github.com/Joseph 
bakulikira/Traveling-Salesman- 
Algorithm 

(Accessed on 21st July 
2023). 


Table 1: Control Variables 
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3.1.2 Experiment 


The following steps were used to obtain a set of results: 


1. 


Run the Brute Force Algorithm with 4 data points and record the shortest distance 
calculated (see Appendix 2). 

Run the elitist GA using the same 4 data points and record the shortest distance 
calculated and the time taken (see Appendix 3). 

Repeat step 2 using ACS (see Appendix 4). 

Repeat steps 2 and 3 two more times to collect repeat readings. 

Repeat all steps using a larger data set (see Appendix 5). 


Calculate the mean time taken for each set of results. 


The default timer() function was imported from the timeit module in Python to obtain the 


time taken for the optimal solution to be calculated as seen in Figure 7. 


t: default timer timer 
Figure 7: Code to import default timer() 
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After program execution, the time and shortest distance were outputted — see Figures 8 
and 9 respectively. 


Time in seconds: 1.5134142320000006 


Figure 8: Example time output 


Genetic Algorithm 


Best distance : 543.89 


Figure 9: Example best distance output 
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3.2 Results 


Measured using the timer function in Figure 7, the results in Tables 2 and 3 show the time 


taken to calculate the shortest distance for each algorithm and data set. The mean time 


taken was found using Equation (6): 


time 1 +time 2+time 3 
3 


mean time = 


An example calculation using Equation (6) is shown below. 


0.233+0.214+0.196 
3 


= 0.2143 


= 0.214 (3 s. f) 


Number of Time taken to find the shortest distance (s) Mean 
Cities Attempt 1 Attempt 2 Attempt 3 time (s) 

4 0.233 0.214 0.196 0.214 

8 1.554 1.219 0.805 1.193 

12 12.02 11.68 12.76 12.15 

16 21.43 19.33 22.46 21.07 

20 29.04 30.21 28.65 29.30 


Table 2: Time taken to find shortest distance using the Elitist GA 


Number of Time taken to find the shortest distance (s) Mean 
Cities Attempt 1 Attempt 2 Attempt 3 time (s) 

4 0.246 0.178 0.197 0.207 

8 0.565 0.689 0.713 0.656 

12 4.595 5.103 4.992 4.897 

16 7.456 7.685 8.103 7.748 

20 12.31 13.85 13.59 13.25 


Table 3: Time taken to find shortest distance using ACS 
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The results from Tables 2 and 3 were plotted to produce Figure 10. 


A graph to show the effect of the number of cities in the tour on the time taken to obtain the shortest 
route 
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Figure 10: A graph to show mean time taken (s) against number of cities 


To evaluate the accuracy of each algorithm, the shortest distance obtained by each of the 


algorithms was noted down in Tables 4 and 5. 


Number of Shortest distance obtained Mean shortest 
Cities Attempt 1 Attempt 2 Attempt 3 distance 
4 543.890 543.890 543.890 543.890 
8 1121.46 1121.46 1121.46 1121.46 
12 1253.88 1261.41 1247.21 1254.17 
16 1531.68 1542.35 1546.54 1540.19 
20 2035.21 1986.42 1964.53 1995.39 


Table 4: Shortest distance obtained using the Elitist GA 


Number of Shortest distance obtained Mean shortest 
Cities Attempt 1 Attempt 2 Attempt 3 distance 
4 543.890 543.890 543.890 543.890 
8 1121.46 1121.46 1121.46 1121.46 
12 1242.87 1239.34 1242.87 1241.69 
16 1515.34 1513.43 1511.76 1513.51 
20 1687.48 1693.78 1701.54 1694.27 


Table 5: Shortest distance obtained using ACS 
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The actual shortest distance was found using the control Brute Force Algorithm (see 


Appendix 2). 
Number of Cities Actual shortest distance 
4 543.890 
8 1121.46 
12 1237.83 
16 1493.76 
20 1632.06 


Table 6: Actual shortest distance obtained using Brute Force 


The data in Tables 4-6 were processed using Equation (5), enabling the percentage 
accuracy of each algorithm to be found (see Table 7). The example calculation shows 


how the percentage accuracy of ACS was determined for 16 cities: 


1493.76 
1513.51 


x 100 = 98.69508626 


= 98.70% (4 s. f) 


Number of Cities Elitist GA Percentage ACS Percentage 
Accuracy (%) Accuracy (%) 
4 100.0 100.0 
8 100.0 100.0 
12 98.70 99.69 
16 96.99 98.70 
20 81.79 96.33 


Table 7: Percentage accuracy calculated using Equation (6) 


The data from Table 7 were plotted to produce Figure 11. 


A graph to show the effect of the number of cities in the tour on the percentage accuracy of each 
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Figure 11: A graph to show percentage accuracy against number of cities 


3.3 Interpretation 


Figure 10 shows that for both algorithms, the time taken to find the optimal route increased 
exponentially as the size of the city tour increased. However, the time taken by the Elitist 
GA increased at a faster rate than ACS suggesting it generally has a greater time 
complexity. When the algorithms navigated through a map of 4 cities or less, both 
algorithms had a similar time complexity as it roughly took them a mean time of 0.2 


seconds for a solution to be generated. 


There were no distinct anomalous results in Tables 2 and 3 since the points in Figure 10 
were all located close to the line of best fit. Extrapolation was used to create a curve 


between 0 and 4 cities for both algorithms, therefore, there was a chance that these points 
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could be unreliable. Nevertheless, this was most likely not the case since background 


theory supports the claim that the time complexity will increase with city tour size. 


Figure 11 clearly illustrates that as the number of cities in the tour increased, the accuracy 
of the distance generated decreased. Both algorithms are proven to be 100% accurate 
for tours that consisted of 8 cities or less; if the city tour was greater, the accuracy 
decreased at a non-linear rate. When experimenting with 20 cities, the Elitist GA had an 
accuracy rate of 81.79% which is significantly lower than the accuracy of ACS which was 
96.33%. Overall, the Elitist GA was shown to be more inaccurate in comparison to ACS 
as the percentage accuracy fell by a faster rate, especially when the city tour was greater 


than 16. 


It is also evident that there were no anomalous results in Tables 4 and 5 as the points in 
Figure 11 only deviated slightly from those on the line of best fit, hence, the results were 
valid and reliable overall. However, the values obtained from the Elitist GA are located 
further away from its line of best fit in contrast to ACS. This implies that the results from 
Table 4 had a greater error uncertainty than those in Table 5. Despite this, the error 
uncertainty did not have an effect on the overall relationship between the accuracy of 


each algorithm and the number of cities. 
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4. Conclusion 


4.1 Research Question Analysis 


This exploration aimed to utilise the theory behind metaheuristic algorithms to formulate 
an experiment, applicable to evaluate each algorithm's time complexity and percentage 
accuracy when solving an NP-hard problem. This investigation showed that ACS has a 
lower time complexity but a higher level of accuracy in comparison to the Elitist GA when 


finding the shortest route to solve the TSP. 


It was proven that accuracy is inversely proportional to the city tour size since larger data 
sets lead to a lower percentage accuracy. On the other hand, the time taken for either 
algorithm to obtain a solution increased with the number of cities visited in the tour. A 
greater time complexity represents a lower level of efficiency; therefore, reinforcing the 


inversely proportional nature between efficiency and data set size. 
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4.2 Hypothesis Analysis 


The hypothesis made was supported by the results to a partial extent as ACS was shown 
to be more accurate, however, the Elitist GA actually has a higher time complexity than 
ACS. After further research, relevant supporting evidence was found. After every iteration, 
the entire population (city tour) changes with ACS whereas with the Elitist GA, only one 


city is altered enabling the optimal solution to be found quicker. 38 


This paper's results are reproducible as a similar experiment, carried out by Alexander, 
A. and Sriwindono, H. had the same findings.? As shown in Figure 12, the distance 
obtained by ACS is always much shorter than that found by the GA, implying that ACS is 
more accurate. The accuracy of the GA tended to decrease at a faster rate than ACS like 
in this paper. 

mACS mGA 


180000 
160000 
140000 
120000 
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40000 all i 
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0 cd A 
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Kota Kota Kota Kota Kota Kota Kota Kota Kota Kota 
mACS 10982 14933 20025 25348 29131 34963 33424 37512 45054 41183 
mGA 13477 29393 46509 64355 83205 100648 112482 132598 147401 167642 


Figure 12: Distance found by each of the two algorithms with shortest distance 
on the y-axis and city tour size on the x-axis?’ 


36 Alexander, A. and Sriwindono, H. (2019) ‘The comparison of genetic algorithm and ant colony optimization in 
completing travelling salesman problem’, Proceedings of the 2nd International Conference of Science and 
Technology for the Internet of Things, ICSTI 2019, September 20th, Yogyakarta, Indonesia [Preprint]. 
doi:10.4108/eai.20-9-2019.2292121. 

37 Ibid. 

38 Ibid. 
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Alexander, A. and Sriwindono, H also recorded the time taken by each algorithm — this is 


shown in Figure 13. 


mACS mGA 


à atl 


10 Kota 20 Kota 30 Kota 40 Kota 50 Kota 60 Kota 70 Kota 80 Kota 90 Kota 


Kota 
mACS 395 125.8 
mGA 4 140.8 


Figure 13: Average time taken by each of the two algorithms with time, in 
seconds, on the y-axis and city tour size on the x-axis *9 


Their data supported the concept that ACS is faster than the GA when navigating through 
a city tour consisting of 20 or fewer. However, when the city tour size exceeds 20, the GA 
was faster than ACS. Therefore, to improve this investigation, data sets larger than 20 


should be used to see whether ACS will have a greater time complexity than GA. 


3? Alexander, A. and Sriwindono, H. (2019), op. cit. 
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4.3 Relevance of Data 


Both algorithms are powerful metaheuristic methods of solving NP-hard problems, so 
comparing, and evaluating the algorithms are beneficial to see which one would 
outperform the other. This experiment proved that ACS is superior to the Elitist GA when 
solving the TSP, a notable optimisation problem. Based on the results and analysis 
presented in this paper, for cities consisting of 20 or fewer, ACS always produced a 


shorter distance in a quicker time period. 


These metaheuristic approaches are applicable beyond the TSP as they can be used to 
solve other NP-hard problems including the Knapsack Problem, Job Shop Scheduling 
Problem and Vehicle Routing Problem. They can approximate a solution more efficiently 


than exact algorithms; hence, this data is insightful. 


Outside of computational complexity theory, these algorithms have real-life applications, 
for instance, in telecommunications network optimisation and circuit board manufacturing. 
In telecommunications, metaheuristic algorithms are used to determine how data packets 
should be routed to improve network performance and reduce latency. On the other hand, 
when manufacturing circuitry, the optimal component arrangement is vital to eliminate 
signal interference. For these reasons, this data is relevant as it can be applied to these 


scenarios to save resources, such as money, and increase efficiency. 
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4.4 Evaluation 


Seeing that the findings of this experiment aligned with the results of other research 


papers (see Section 4.2), this essay was proven to be successful and had many 


strengths. The accuracy of the results was ensured through repeat readings and 


calculating a mean. In addition, the absence of anomalies supports emphasises the 


validity of the results. Despite this, several improvements would eliminate potential 


sources of error since the investigation was carried out in a home environment — these 


are presented in Table 8. 


Source of Errors 


Effect and Importance 


Improvements 


Brute Force Algorithm: For 
large data sets, the time 


taken to generate a solution 
was extremely slow since it 
has a time complexity of 
O(n!). 


The program may have 
stopped too soon, meaning 
the final solution generated 
was inaccurate. As this 
value was used to calculate 
the percentage accuracy, 
this would have had an 
effect on the values in 
Table 15. 


Other exact algorithms with 

a lower time complexity 

could be used: 

- Brand and 
0(2”) 

- Dynamic Programming: 
O(n? x 2”) 


Bound: 


Random Time Error: The 
programs’ performance can 
be affected by other 
computer processes 
running simultaneously. 


Although all other 
applications were closed 
during program execution, 
background processes are 
still being carried out. This 
would affect the accuracy 
of the default timer(), 
impacting the results in 
Tables 2 and 3. 


A computer with a better 
processor and larger RAM 
could eliminate this source 
of error. 


Systematic Error from 
Source Code: The program 


code was found online so 
there may be logic errors. 


This error would have an 
effect on most of the results 
obtained if the code does 
not directly mirror the 
algorithms. 


One could program their 
own algorithms or use code 
that they are certain to 
contain no logic errors. 


Table 8: Potential errors and corresponding improvements 
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After evaluating the data acquired in this investigation, it has been concluded that ACS is 
both more accurate and has a lower time complexity than the Elitist GA when it comes to 
solving the TSP for 20 cities or fewer. Therefore, the research question “How does the 
Elitist Genetic Algorithm compare to Ant Colony System in terms of time 
complexity and accuracy when attempting to solve the Travelling Salesman 
Problem?” has been answered, consequently underscoring the beneficial nature and 


indispensability of this essay. 
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6. Appendix 


All the code was accessed on 21% July 2023 from https://github.com/Joseph 
bakulikira/Traveling-Salesman-Algorithm. See the following pages for screenshots of the 


code used. 
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6.1 Appendix 1 — Main Code 


main.py 
import 

from import * 

from import * 

from import 

from : import * 

from import SumDistance 


+ 


from import defa 
.init() 


0 
es = ["ACS", "ELITIST", "MAX-MIN"] 


manager. Background() 


delta time = anager.SetFps() 
manager. UpdateCaption() 


.get(): 


.QUIT: 


j r() 
print('Time in seconds: 
if event.key == .K SPA 
| " 


e .MOUSEBUTTONDOWN : 
button -- 1: 


htMouseClicked 


manager.DrawPoints() 
manager.DrawShortestPath() 


elif selectedIndex == 
if pause == E 
manager.Lexicographic() 
manager.DrawPoints() 
manager .DrawShortestPath() 
manager.Percentage(manager.PossibleCombinations) 
elif selectedIndex -- 
if pause == E 
manager.GeneticAlgorithm() 
manager.DrawPoints() 
manager.DrawShortestPath() 
else: 
manager.AntColonyOptimization(pause) 


manager.ChangeAntColonyVariation(antColonyTypes [selectedIndex-3]) 
manager.Percentage((literations)) 
manager.ShowText(selectedIndex, started) 


if showUI: 
panel.Render(manager.screen) 
AlgorithmChoice.Render(manager.screen, rightMouseClicked) 
if pause !- PauseButton.state: 
PauseButton.state - pause 


PauseButton.Render(manager.screen, rightMouseClicked) 
ResetButton.Render(manager.screen, rightMouseClicked) 
RandomButton.Render(manager.screen, rightMouseClicked) 


pause PauseButton.state 
reset ResetButton.state 


if reset == 
reset = 
ResetButton.state = 
temp = manager.Points.copy() 
manager = (temp) 
manager.OptimalRoutes = manager.Points.copy() 
manager.recordDistance = SumDistance(manager.Points) 
manager.ResetAntColony(manager.antColony.variation) 
manager.ResetGenetic() 


GenerateToggle - RandomButton.state 
if GenerateToggle == 
manager.RandomPoints () 
GenerateToggle - 
RandomButton.state 
if pause -- - 
PauseButton.text = "Continue" 
else: 
PauseButton.text = "Pause" 
if rightMouseClicked: 
selectedIndex = AlgorithmChoice.currentIndex 


manager.scaler += 1 
if manager.scaler > manager.max radius: 
manager.scaler = manager.max radius 
.flip() 
rightMouseClicked = 
.quit() 


6.2 Appendix 2 — Brute Force Algorithm Code 


manager.py > Y Manager > ® ResetAntColc 


import 
import 


from 
from 
from 
from 
from 


from 


import randint 
import 
import * 
import 
import x 
import 


100 
= 1920, 1080 
= 150 


= 150 
= 300 
«init() 


(width, height) 
30 
»Set mode(size) 
0 


algorithms = ["Brute Force", "Lexicographic Order", "Genetic Algorithm", "Ant Colony ACS", "Ant Colony Elitist", "Ant Colony Max-Min"] 


genetic ( [sample ( (n)), n) for populationSize)], populationSi 


PossibleCombinations = Factorial(n_points) 
print("possible combinations : ".format(Factorialí(n points 


Order i for (n points)] 


— init (self): 
Points [ (372, 201), (298, 247), (187, 189), (458, 325), (398, 289), (198, 456), (346, 150), (162, 157), 
self.Points = Points 
self.reco SumDistance(self.Points) 
1 Points.copy() 
oints.copy() 


self.antColony - (variation-"ACS", | size, max ite 


nodes=self.Points.copy(), a beta rho=0.1, pheromon 


ResetGenetic(self): 


self.genetic = ( [sample (n)), n p ionSize)], populationS 


ResetAntColony(self, name="ACS"): 

self.recordDistance = SumDistance|[self.Points|] 

self.antColony - (variation-name, size-colony size, max iterations - iterations, 
nodes=self.Points.copy(), alpha=1, beta=3, rho=0.1, pheromone-1, phe deposit weight-1) 

SetFps(self): 

return self.clock.tick(self.fps)/1000.0 


UpdateCaption(self): 
frameRate - (self.clock.get fps()) 
à .set caption("Traveling Salesman Problem - Fps : " , format(frameRate)) 


Counter(self): 

self.counter += 1 

if self.counter » self.PossibleCombinations: 
self.counter = self.PossibleCombinations 


BruteForce(self): 
if self.counter != self.PossibleCombinations: 
il randint(0, self.n_points-1) 
i2 randint(0, self.n points-1) 
self.Points[i1], self.Points[i2] = self.Points[i2], self.Points[i1] 


(201, 317)] 


dist = SumDistance(self.Points) 

if dist < self.recordDistance: 
self.recordDistance = dist 
self.OptimalRoutes = self.Points.copy() 


self.DrawLines() 

Lexicographic(self): 

self.Order = LexicalOrder(self.Order) 

nodes = [] 

for i in self.Order: 
nodes.append(self.Points[il) 


self.Counter() 


dist = SumDistance(nodes) 

if dist « self.recordDistance: 
self.recordDistance - dist 
self.OptimalRoutes = nodes.copy() 


self.DrawLines() 


GeneticAlgorithm(self) 
self.genetic.CalculateFitness(self.Points) 
self.genetic.NaturalSelection() 


for i in (self.n_points): 
self.currentList[i] = self.Points[self.genetic.current[i]] 
if self.genetic. record < self.recordDistance: 
for i in (self.n_points): 
self .OptimalRoutes[il = self.Points[self.genetic.fitest[i]] 
self. recordDistance = self.genetic. record 


self .DrawLines( 


AntColonyOptimization(self, pause): 
if pause == : 
self.counter 4- 1 
if self.counter » self.antColony.max iterations: 
self.counter = self.antColony.max iterations 


if self.counter « self.antColony.max iterations: 
self.antColony.Simulate(self.counter) 


self.antColony.Draw(self) 
self.recordDistance - self.antColony.best distance 


RandomPoints(self): 

self.Points - [[ (167,92), (73,5), (56,24), (6,41) 
self.ResetAntColony(self.antColony.variation) 

self.recordDistance = SumDistance(self.Points) 

self.OptimalRoutes = self.Points.copy() 

self.currentList = self.Points.copy() 


Percentage(self, val): 
percent = (self.counter/val) x 100 
textColor = (255, 255, 255) 


textFont = ; .SysFont("Arial", 20) 
textSurface = textFont.render( (round(percent, 4)), , textColor) 
self.screen.blit(textSurface, (width//2, 50)) 


def ShowText(self, selectedIndex, started = True): 
textColor = (255, 255, 255 


textFont = pygame.font.SysFont 
textFont2 = pygame.fo SysFont 


textSurfacel - textFont.render * str(round(self.recordDistance,2)), False, textColor 
textSurface2 = textFont.render(self.algorithms [selectedIndex], False, textColor 
textSurface3 textFont2.render 3 ‚False, textColor 


self.screen.blit(textSurface1, (100, 70) 
self.screen.blit(textSurface2, (100, 35) 
if started == False: 

self.screen.blit(textSurface3, (width//2, height-200) 


def DrawShortestPath(self): 
if len(self.OptimalRoutes) > 0: 
for n in range(self.n_points): 
_i = (n+1)%self.n_points 
pygame.draw.line(self.screen, self.Highlight, 
(self .OptimalRoutes [n].x, self.OptimalRoutes[n].y), 


(self .OptimalRoutes[_il.x, self.OptimalRoutes[ il.y), 
self.LineThickness 
self .OptimalRoutes[n].Draw(self, self.showIndex, True, n 


def DrawPoints(self, selected_index = 0): 
for point in self.Points: 
point.radius = self.scaler 
point.Draw(self 


def DrawLines(self, drawCurrent=False 
if drawCurrent == True: 
for i, point in enumerate(self.currentList): 
_i = (i+1)%len(self.currentList 
pygame.draw.line(self.screen, self.Gray, (point.x, point.y), (self.currentList[ il.x, self.currentList[ il 
else: 
for i, point in enumerate(self.Points): 
_i = (i+1)%len(self.Points 
pygame.draw.line(self.screen, self.Gray, (point.x, point.y), (self.Points[_il.x, self.Points[_il.y), 1 


def Background(self): 
self .screen.fill(self.Black 
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6.3 Appendix 3 — Elitist Genetic Algorithm Code 


genetic_algorithm_TSP.py > ® createRoute 
import as A 5 A 


init (self, xy). 
self.x 2 x 
self.y = y 


distance(self, city): 

xDis = abs(self.x - city.x) 

yDis = abs(self.y - city.y) 

distance = np.sqrt((xDis ** 2) + (yDis xx 2)) 
return distance 


_repr_ (self): 
return "(" + (self.x) + "," + (self.y) + ")" 


. init (self, route): 
self. route = route 
self .distance = 0 
self.fitness= 0.0 


routeDistance(self): 
if self.distance ==0: 
pathDistance = 0 
for i in (0, len(self.route)): 
fromCity = self.routeli] 
toCity = 
if i + 1 < len(self. route): 
toCity = self.route[i + 1] 
else: 
toCity = self.route[0] 
pathDistance += fromCity.distance(toCity) 
self.distance - pathDistance 
return self.distance 


routeFitness(self): 
if self.fitness -- 0: 

self.fitness = 1 / (self.routeDistance()) 
return self.fitness 


createRoute(cityList): 
route - .sample(cityList, len(cityList)) 
return route 


initialPopulation(popSize, cityList): 
population - [] 


for i in (0, popSize): 
population.append(createRoute(cityList)) 
return population 


rankRoutes (population): 
fitnessResults - () 
for i in (0, len(population)): 
fitnessResults[i] = (population[il).routeFitness() 
return sorted(fitnessResults.items(), key - : (1), reverse 


selection(popRanked, eliteSize): 

selectionResults [] 

df = pd. (np.array(popRanked), columns=["Index","Fitness"]) 
df['cum sum'] = df.Fitness.cumsum() 

df['cum perc'] = 100xdf.cum sum/df.Fitness.sum() 


for i in (0, eliteSize): 
selectionResults.append(popRanked[il [0]) 
for i in (0, len(popRanked) - eliteSize): 
pick = 100x . random( ) 
for i in (0, len(popRanked)): 
if pick <= df.iat[i,3]: 
selectionResults.append(popRanked [i] [0]) 
break 
return selectionResults 


matingPool(population, selectionResults) : 

matingpool = [] 

for i in (0, len(selectionResults)): 
index = selectionResults[i] 
matingpool.append(population[index]) 

return matingpool 


breed(parent1, parent2): 
child = [] 

childP1 - [] 

childP2 = [] 


geneA .random() * len(parent1)) 
geneB .random() x len(parent1)) 


startGene = min(geneA, geneB) 
endGene = max(geneA, geneB) 


for i in (startGene, endGene): 
childP1.append(parent1[i]) 


childP2 = [item for item in parent2 if item childP1] 


child = childP1 + childP2 
return child 


breedPopulation(matingpool, eliteSize): 

children = [] 

length = len(matingpool) - eliteSize 

pool - .sample(matingpool, len(matingpool)) 


for i in (0, eliteSize): 
children.append(matingpool[il) 


for i in (0, length): 
child = breed(pool[il, pool[len(matingpool)-i-1]) 
children.append(child) 

return children 


mutate(individual, mutationRate): 
for swapped in (len(individual)): 
if ( .random() « mutationRate) : 
swapWith = ( .random() * len(individual)) 


cityl = individual [swapped] 
city2 = individual [swapwWith] 


individual[swapped] = city2 
individual [swapWithl = city1 
return individual 


mutatePopulation(population, mutationRate): 
mutatedPop - [] 


for ind in (0, len(population)): 
mutatedInd - mutate(population[ind], mutationRate) 
mutatedPop.append(mutatedInd) 

return mutatedPop 


nextGeneration(currentGen, eliteSize, mutationRate): 
popRanked - rankRoutes(currentGen) 

selectionResults - selection(popRanked, eliteSize) 
matingpool = matingPool(currentGen, selectionResults) 
children = breedPopulation(matingpool, eliteSize) 
nextGeneration = mutatePopulation(children, mutationRate) 
return nextGeneration 


geneticAlgorithm(population, popSize, eliteSize, mutationRate, generations): 
pop = initialPopulation(popSize, population) 
print("Initial distance: " + (1 / rankRoutes(pop) [0] [11)) 


for i in (0, generations): 
pop = nextGeneration(pop, eliteSize, mutationRate) 


print("Final distance: " + (1 / rankRoutes(pop) [0] [11)) 
bestRouteIndex = rankRoutes (pop) [0] [0] 

bestRoute = pop[bestRouteIndex] 

return bestRoute 


cityList = [] 


for i in (0,25): 
cityList.append( (x= .random() * 200), y= .random() * 200))) 


print(cityList) 


geneticAlgorithm(population=cityList, popSize=150, eliteSize=20, mutationRate=0.01, generations=500) 


geneticAlgorithmPlot(population, popSize, eliteSize, mutationRate, generations): 
pop = initialPopulation(popSize, population) 

progress = [] 

progress.append(1 / rankRoutes(pop) [0] [1] ) 


for i in (0, generations): 
pop = nextGeneration(pop, eliteSize, mutationRate) 
progress.append(1 / rankRoutes(pop) [0] [11) 


.plot(progress) 
.ylabel('Distance') 
.Xlabel('Generation') 
. show( ) 
geneticAlgorithmPlot(population-cityList, popSize-100, eliteSize-20, mutationRate-0.01, generations-500) 


from utils import x 
from random import randint 


class Genetic: 

def ^ init (self, population=[], populationSize-0): 
self.population - population 
self.size - populationSize 
self.fitness = [0 for i in range(populationSize) 
self.record - float 
self.currentDist - float 
self.current - None 
self.fitest = 
self.fitestIndex - 0 
self.mutation rate - 0.01 


CalculateFitness(self, points): 
for i in range(self.size): 
nodes = 
for j in self.population[i 
nodes.append(points[j] 


dist = SumDistance(nodes 
if dist « self.currentDist: 
self.current = self.population[i 


if dist « self.record : 
self.record - dist 
self.fitest = self.population[i 
self.fitestIndex - 


self.fitness[i] = 1/ (dist+1 
self.NormalizeFitnesss 


NormalizeFitnesss(self): 

s=0 

for i in range(self.size): 
S += self.fitness[i 

for i in range(self.size): 
self.fitness[i] = self.fitness[il/s 


Mutate(self, genes): 
for i in range(len(self.population[0]) 
if (randint(0, 100)/100) « self.mutation rate: 
a - randint(0, len(genes)-1 
b = randint(0, len(genes)-1 
genes[a], genes[b] = genes[b], genes [a 


CrossOver(self, genesi, genes2): 
start = randint(®, len(genes1)-1 
end = randint(start-1, len(genes2)-1 
try: 
end = randint(start+1, len(genes2)-1 
except: 
pass 
new_genes = genes1[start:end 
for i in range(len(genes2) 
p = genes2li 
if p not in new_genes: 
new_genes.append(p 


return new_genes 


NaturalSelection(self): 

nextPopulation = 

for i in range(self.size): 
generation1 = PickSelection(self.population, self.fitness 
generation2 = PickSelection(self.population, self. fitness 
genes = self.CrossOver(generationi, generation2 
self.Mutate(genes 
nextPopulation.append(genes 

self.population = nextPopulation 
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antColony.py > ... 


import 

from import * 

from import translateValue 
o .init() 

textColor - (0, 0, 0) 


textFont .SysFont("Arial", 20) 


( Je 
. init (self, variation-"ACS", size=5, elitist weight-1.0, minFactor-0.001, alpha=1.0, beta=3.0, 
rho=0.1, phe deposit weight-1.0, pheromone=1.0, max iterations-100, nodes- ; E 
„variation = variation 
„size = size 
.elitist weight = elitist weight 
.minFactor = minFactor 
.alpha = alpha 
.rho = rho 
.phe deposit weight - phe deposit weight 
f.max iterations - max iterations 
f.n nodes = len(nodes) 
f.nodes = nodes 
.edges = [[ for j in self.n_nodes)] for i in (self.n_nodes) ] 
for x in (self.n_nodes): 
for y in (self.n_nodes): 
heuristic = .sqrt( 
.pow(self.nodes [x].x-self.nodes [y].x, 2) + 
.pow(self.nodes [x].y-self.nodes [y].y, 2) 


) 
self.edges[xl[y] = self.edges [y] [x] = (x, y, heuristic, pheromone) 
feants = [ (self.edges, alpha, beta, self.n_nodes) for i in (self.size)] 


lf.best tour = [] 
f.best_distance = ("inf") 


self.local best route = [] 
self.local best distance = ("inf") 


AddPheromone(self, tour, distance, heuristic-1) 
= self.phe deposit weight / distance 
for i in (self.n nodes): 
self.edges[tour[ill[tour[(i + 1) % self.n nodes]l.pheromone += heuristic 


ACS(self): 


for ant in self.ants: 
self.AddPheromone(ant.UpdateTour(), ant.CalculateDistance()) 
if ant.distance « self.best distance: 
self.best tour - ant.tour 
self.best distance - ant.distance 


for x in (self.n nodes): 
for y in (x * 1, self.n nodes): 
self.edges[x] [y]. pheromone *= (1 - self.rho) 


Simulate(self, counter): 

if self.variation == "ACS": 
self.ACS() 

elif self.variation == "ELITIST": 
self.ELITIST() 

elif self.variation == "MAX-MIN": 
self.MAX MIN(counter) 


Draw(self, manager): 


for i (len(self.best tour)): 

self.nodes[self.best tour[i]] 
self.nodes[self.best_tour[(i+1) % len(self.best_tour)]] 
.line(manager.screen, manager.Highlight, a.GetTuple(), b.GetTuple(), manager.LineThickness) 


variation == "MAX-MIN": 
ant in self.ants: 
for edge in ant.edges: 
for e in edge: 
b (min((e.pheromone)x: .pow(10, 5), 255)) 
thickness (translateValue(e.pheromone, 0, 255, 1, 8)) 


.line(manager.screen, (r, g, 0), self.nodes[e.a].GetTuple(), self.nodes[e.b].GetTuple(), thickness) 


in self.ants: 

edge in ant.edges: 

for e in ed 
(min((e.pheromone)*2, 255)) 
(translateValue(e.pheromone, 0, 255, 1, 8)) 


.line(manager.screen, (r, g, 0), self.nodes[e.a].GetTuple(), self.nodes[e.b].GetTuple(), thickness) 


for node in self.nodes: 
.circle(manager.screen, manager.White, node.GetTuple(), manager.scaler) 


for i in self.best tour: 


textSurfa textFont.render( (i), , textColor 


textRectangle = textSurface.get_rect(center=(self.nodes[self.best_tourli]].x, self. 


ger.screen.blit(textSurface, textRectangle) 


ant.py > ... 
import 
import 
import 


__init__(self, a, b, heuristic, pheromone): 
self.a=a 

self.b = b 

self.heuristic heuristic 

self.pheromone = pheromone 


__init_ (self, edges, alpha, beta, n nodes): 


self.edges edges 
self.tour 

self.alpha - alpha 
self.beta - beta 
self.n nodes - n nodes 
self.distance - 0.0 


nodes [self. 


t tour[i] 


.y)) 


def NodeSelection(self): 
Constructing solution 
an ant will often follow the st 
ne trail v 
point on raph or 
ecting the next city depending on the d 


e amount of pheromone on the path be 


"i 


roulette wheel - 0 
states - [node node in range(self.n nodes) if node not in self.tour] 
heuristic value = 0 
new state in states: 
heuristic value += self.edges[self.tour[-1]] [new state].heuristic 
new state states: 
A = math.pow(self.edges[self.tour[-1]] [new statel.pheromone, self.alpha) 
B = math.pow((heuristic value/self.edges[self.tour[-1]] [new statel.heuristic), self.beta) 
roulette wheel += A x B 
random value - random.uniform(0, roulette wheel) 
wheel position - 0 
new state states: 
A = math.pow(self.edges[self.tour[-1]] [new statel.pheromone, self.alpha) 
B = math.pow((heuristic value/self.edges[self.tour[-1]] [new statel.heuristic), self.beta) 
wheel position += A * B 
wheel position »- random value: 
new state 


def UpdateTour(self): 
self.tour - [random.randint(0, self.n nodes - 1)] 
len(self.tour) « self.n nodes: 
self.tour.append(self.NodeSelection()) 
self.tour 


def UpdateTour(self): 
self.tour - [random.randint(0, self.n nodes - 1)] 
len(self.tour) « self.n nodes: 
self.tour.append(self.NodeSelection()) 
self.tour 


def CalculateDistance(self): 
self.distance - 0 
i in range(self.n nodes): 
self.distance += self.edges[self.tour[i]] [self.tour[(i+1)%self.n_nodes]].heuristic 
|! self.distance 
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6.5 Appendix 5 — Data Sets 


6.5.1 Four City Map 
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6.5.3 Twelve City Map 
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6.5.4 Sixteen City Map 
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6.5.4 Twenty City Map 
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